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Abstract 

We consider the most general loop integral that appears in non-relativistic effective field 
theories with no light particles. The divergences of this integral are in correspondence with 
simple poles in the space of complex space-time dimensions. Integrals related to the original 
integral by subtraction of one or more poles in dimensions other than D = 4 lead to nonminimal 
subtraction schemes. Subtraction of all poles in correspondence with ultraviolet divergences 
of the loop integral leads naturally to a regularization scheme which is precisely equivalent to 
cutoff regularization. We therefore recover cutoff regularization from dimensional regularization 
with a nonminimal subtraction scheme. We then discuss the power-counting for non-relativistic 
effective field theories which arises in these alternative schemes. 
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1 Divergent integrals in non-relativistic effective field theories 



Effective field theories of non-relativistic scattering have been a subject of much interest recently. 
These field theories have applications in many branches of physics. They are useful in analyzing 
a system with widely separated energy scales. For instance, to describe two-body scattering at 
momenta considerably lower than the mass of any exchanged particle an effective field theory can 
be written in which all of the exchanged particles are "integrated out" of the original field theory 
Lagrangian. In the low-energy effective theory the effects of these "heavy" degrees of freedom 
are represented by an expansion of the interaction Lagrangian as a sequence of local operators of 
increasing dimension. This low-energy effective theory can then be used to calculate non-relativistic 
scattering in the system of interest. 

In pursuing such a calculation one encounters divergences. These divergences must be regulated 
and renormalized before physical quantities can be calculated. Here we consider the most general 
loop integral appearing in a low-energy EFT description of non-relativistic scattering. We regularize 
this integral using the standard technique of dimensional regularization (DR). Within the framework 
of DR, we identify the poles as a function of the number of space-time dimensions that correspond to 
the divergences of this integral. This allows us to clarify the relation between DR-based schemes and 
cut-off regularization. We also demonstrate explicitly how the power-counting for the coefficients 
appearing in the non-relativistic effective theory Lagrangian is affected by the choice of regulator. 

Consider a non-relativistic effective field theory for the interaction of two identical heavy par- 
ticles of mass M. If no exchanged degrees of freedom appear explicitly in the Lagrangian then we 
can immediately write: 

C = ^id t ip + - \co{^H? - \c 2 {^V 2 ^H) + h.c. + .... (1) 

where the dots refer to other invariants with two or more derivatives. Of course, without a power- 
counting scheme this Lagrangian is useless, since there are infinitely many interaction terms con- 
sistent with assumed symmetries. A power-counting scheme must be established which allows us 
to keep only a finite number of operators at a given order in the gradient expansion. Unless this 
is done the effective field theory will have no predictive power, as an infinite number of coefficients 
will enter the calculation. 

We can formally write the all-orders solution of the low-energy effective theory by constructing 
the non-relativistic potential 

oo 

V(p) = £ C 2n p 2n (2) 

n=0 

directly from the Lagrangian. Here p is to be understood as an operator, which may denote 
momentum or energy-dependence. More explicitly, this potential is 

oo n 

F(k',k;£) = EEE*V M ^(k',k), (3) 

n=0 i=l j 

where p 2 = ME and the index j enumerates members of the set of all operators of dimension 2i 
which are consistent with Hermiticity and rotational invariance. In the case of s-wave scattering, 
the operators that contribute are 

O) = k 2j k' 2(i ~ j) + k 2ii - j) k' 2j , (4) 
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although, as we shall see, these details of the construction of the potential are not necessary for our 
discussions here. Note that the potential V contains only terms that are analytic in the momenta 
k, k' and p. This is because in writing the effective Lagrangian (|l|) we have made local expansions 
of any non-analytic structures which appear in the "full" theory. Therefore, the only non-analytic 
effects in this low-energy effective theory come from solving the Schrodinger equation. The effective 
theory is therefore only valid for energies well below the masses or production thresholds of any 
exchanged particles. 

The potential V is then iterated via the Lippmann-Schwinger equation 

T(k\ k; E) = V(k', k; E) + M J V(k', q; £) — - . T(q, k; E) (5) 

to give the on-shell scattering amplitude T(p,p; E), with p = V ME. All observables in the tptp 
system can be obtained from this amplitude. In the language of Feynman diagrams, the Lippmann- 
Schwinger equation generates the sum of all graphs, as illustrated in Fig. |l|. 



Figure 1: The diagrammatic solution of the Lippmann-Schwinger 
equation with the effective potential represented by the shaded 
blob. 

It immediately follows from Eq. (Jsj) that divergent integrals will accompany every loop computed 
in the effective field theory. This means that the C2n coefficients must be renormalized, and 
therefore are renormalization scheme dependent. Since the scaling behavior of the Cm coefficients 
determines which operators are big and which are small in the effective theory, this means that the 
power counting scheme will look different in different regularization schemes. 

In what follows we will consider the most general divergent integral that appears in the effective 
theory. In Section 2 we will explore the many possible definitions of the divergent integral in 
dimensional regularization, and then in Section 3 discuss the power counting for the coefficients 
C2n that arise from these definitions. 

The most general form of the divergent integral which appears in the effective theory is: 

where D is the number of space-time dimensions, \i is a renormalization scale, and E + = E+irj, with 
i] a positive infinitesimal, is the energy. Intuition about this integral may be gained by evaluating 
with a sharp cutoff j3. In D = 4 this gives 

W)=M *** _ a =^/ dq * _ a . (7) 



(2vr) 3 ME+ - q 2 2ir 2 J H ME+ - q 2 ' 

The imaginary part is, of course, regularization scheme invariant. Meanwhile, the real, or principal 
value, part is given by 

Re , M = -iL fa + + ... + ^] + ^o S (^), (8) 



2n+l * Att 2 \/3-p 



provided that (3 2 > p 2 . It is clear that the integral of interest is ultraviolet power-law divergent at 
D = 4. The maximum power of the cutoff which appears is just the superficial degree of divergence 
of the graph, 2n + 1. 

Using the usual techniques employed in calculating dimensionally regularized integrals in an 
arbitrary number of dimensions (see, for instance, Appendix B of Ramond |l|]), the integral (|6|) 
may be re-expressed as, 

, s />\ 4 " D M r ( 2n+i3-i ) /3-D-2n\ 

*M = - (I) (47r p-D/ 2 \x4±T { - ME) / T { — 2 — ) (9) 

provided that E < 0. This equality is valid in a region of the complex D-plane defined by 

1 - 2n < D< 3 - 2n and D > 1. (10) 

The first inequality is easy to understand from Eq. (^) using simple power counting arguments. 
The condition D > 1 — 2n ensures that integral is free of infrared divergences and the condition 
D < 3 — 2n ensures that the integral is free of ultraviolet divergences. According to Eq. (|9|), 
in the complex D-plane (see Fig. ||) these divergences manifest themselves as singularities of the 
gamma functions in the corresponding regions. The condition D > 1 follows from requiring that 
the space-time measure be infrared safe. 




-3-2n 



-l-2n 



5-2n 



l-2n 3-2n 



Figure 2: The complex D-plane. The integral I n {lj) is given 
by Eq. (11) within the shaded region of the D-plane. Poles are 
represented by crosses. 



The two conditions of Eq. ( |10| ) are incompatible for n > 0. Thus, in fact, to properly calculate 
the integral (||) we must do the angular integration in D = 4, and then consider the result of varying 
the number of dimensions in which the "radial" integration is performed. This technique leads to 
the result: 



which is valid for 



4-D 



m r( 



(4vr) 3 / 2 r(§) 



-(-MD) (2n+D_3)/2 r 



3 - D - 2n 



1 - 2n < D < 3 - 2n, 



(11) 



(12) 



and E < 0. If we are interested in scattering problems in four space-time dimensions we must now 
perform two analytic continuations: one in the energy E and one in the number of dimensions D. 
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Performing the analytic continuation in E first, by using the evaluation at E + irj, rather than E 
itself, we find that for E > 



w = (|)'- D ( _M £ ).(-M E+ ,<-3)/^r (*t±^zl) r (*z*z*) , (is, 

provided that 1 — 2n < D < 3 — 2n. The following discussion makes extensive use of this result. 

The usual dimensional regularization prescription for calculation of the integral I n in dimensions 
where the expression (|l~3| ) is not valid is to define a new function I n (/i) which agrees with the original 
I n {n) in the region 1 — 2n < D < 3 — 2n, but is defined by Eq. ( fJ3| j in the rest of the complex 
D-plane. For this I n we may use the Gamma function identity, 

T(n + a)r(l -a-n) = T{a)T{\ - a)(-l) n , (14) 

to obtain 

lM = ^ (f T° (MEn-ME+r-^T (^) r (^) , (15) 

which is correct for all D except odd integers. This expression differs from that of Ref. |2|, |5| by 
a factor due to the way the angular integration was performed. Note that if the identity ( |14| ) is 
applied directly to Eq. then the result of Refs. |2|, [| for I n (fi) is recovered. 



2 Subtraction schemes and the divergent integral 

In D = 4 Eq. (O) yields the result of Ref. [|: 



2n iM P 



InM = ~P 2n ^f , (16) 



where p = y/ME. The regularization scheme corresponding to Eq. (|l~5| ) with removal of any poles 
at the critical dimension is, by convention, minimal subtraction. Of course, there are no poles in 
the dimension of interest to us, D = 4, and so only the imaginary part of the integral survives in 
minimal subtraction. 

The expression ( |T3| ) has poles at D = 1 — 2n, —1 — 2n, . . ., D = 3 — 2n, 5 — 2n, 7 — 2n, 
The former lie below the region of analyticity of I n (see Fig. ^), and can be considered infrared 
divergences that arise when I n is evaluated in such dimensions. Conversely, the latter poles can be 
considered ultraviolet divergences of I n . Some of them lie between the region of analyticity and 
the critical dimension D = 4 which we wish to analytically continue to. In the rest of this work 
we follow a suggestion of Kaplan, Savage and Wise ||, ||] and consider the effects of cancelling out 
poles of r( 3 ~^~ 2n ) (ultraviolet singularities) and r( 2n+ 2 D ~ 1 ) (infrared singularities) in the integral 
I n . This procedure leads to regularization schemes which differ from the standard dimensional 



regularization with minimal subtraction result of Ref. g, Eq. (16). 

Of course, after renormalization the on-shell piece of the amplitude T calculated via Eq. (|5|) 
is necessarily the same in any scheme. Otherwise one of the basic tenets of effective field theory, 
insensitivity to short-distance physics, is violated. So we must ultimately get the same physical 
result regardless of how many poles we subtract 0. This equivalence after renormalization is enforced 

4 This point has been stressed in the context of effective field theories of the nucleon-nucleon interaction by van 
Kolck §. 



by making appropriate choices for the coefficients C2n- Thus, the behavior of the coefficients, and 
hence the power-counting for these coefficients, is affected by the scheme we use to define l n . As 
we will show in the next section some schemes are more useful for reproducing certain physics in 
the effective theory. 

Consider first the gamma function poles corresponding to ultraviolet divergences in cutoff reg- 
ularization. When D is close to 3 — 2m, with m > 0, the poles of r( 3 ~^~ 2n ) between D = 3 — 2n 
and D = 4 have the structure 

t\/t . 2m+l i 
J}LE n 2(n-m) 1 (m 

4tt 2 2 2m 1 D-3 + 2m' { ' 

Thus we can cancel these poles out by defining, 

M n /i,\ 2m+1 1 
/new = j { ) _ M y »\ 2(n-m) 1 

™vp; 2tt 2 ^12/ D-3 + 2m v ; 



m=0 

where I n {^) is given by expression 



Similarly, to cancel the poles of T( - D 2 2n ) in D > 4 dimensions we must add a term 



j^y 2(n + m)(/A 1 (19) 



m=l 



Thus if we wish to cancel all poles of the expression (13) for dimensions D > 3 — 2n we must define 



2vr 2 ^ \2 1 D-3 + 2m' 

m=—oo y ' 

This subtracts the poles which correspond to ultraviolet divergences for D > 3 — 2n. We refer to the 
resulting regularization scheme as dimensional regularization with ultraviolet power-law divergence 
subtraction, or uvPDS. 

It is now a straightforward matter to show that uvPDS is precisely equivalent to cutoff regu- 
larization in D = 4. For the pieces of the sum with m > we find 



M ^ {H\ 2m+ \„ 2(n - m) 1 M_ 



9^2 L^i I o 



2vr 2 ^ n V2y P l + 2m 2vr 2 

m=0 



3 ^2n-2 / ,,\ 2n+l ^ 



|)^ + ( S )'^ + ... + ( f ) 2n + i 

(21) 

and for the pieces which arise from poles in D > 4 it follows from the observation that 

— T P 2m ( ^ m+1 - ^ log (tB±l\ (22) 

2**^* \2) 2m-l-4vr 2i0g U/2-pJ' ^ 

provided that p is kept less than fi/2. It is clear that if we identify j3 = /x/2, then in D = 4 I" ew in 
uvPDS is precisely equivalent to the cutoff result of Eq. (pi). 



3 Scaling of effective field theory coefficients 

The choice of I n made by Kaplan et al, known as power-law divergence subtraction, or PDS, is 
to cancel the pole of expression (15) in D = 3 for all n g, [3|. Of course, as we have seen the 
PDS scheme is only one of a much more general class of definitions of the integral I n . Kaplan et 
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oL's main motivation in considering the PDS scheme was to produce a consistent power counting 
in non-relativistic effective theories where there is an unnaturally large scattering length. Much 
blood, sweat, and ink has been spilt over this problem during the past few years [2-26]. The PDS 
approach to this difficulty is elegant, and leads to good power counting for the coefficients in the 
low-energy Lagrangian. Here we show that in fact PDS is the simplest scheme of this type which 
leads to good power-counting for the bare coefficients C2 n - It has been shown how to obtain similar 
power-counting for the on-shell T- matrix within the framework of any regularization scheme [||, 18] 



However, Refs. (5|, [l8| do not address the issue of the scaling of the coefficients in the bare EFT 
Lagrangian. The power counting of Refs. |2|, [3|, ||, 18] has also been derived within the framework 
of the Wilsonian renormalization group . 

3.1 Review of coefficient scaling with the PDS definition of the divergent inte- 
gral 



The PDS choice is to retain only the m = term of the sum (|20D. This corresponds to the linear 
divergence in the cutoff approach. The integral I n is defined to be (throughout this section we work 
in D = 4): 

I p n^ = -P 2n ^iw+^). (23) 

where our result for the second term differs by a factor of ir from that of Refs. ||, [| because the 
angular integration was performed in D = 4 rather than D = 3. 

This definition of I n leads to a straightforward solution of the Lippmann-Schwinger equation, 

1 1 Ma iMp , 

+ + TT- (24) 



T on (p) J2n=oC 2 nP 2n 4vr 2 4vr 

Now, we attempt to match this to a form of the inverse amplitude which corresponds to the presence 
of an unnaturally large scattering length 

M ( 1 1 9 4n \ 

- + -r e p 2 + 0(p 4 ) - ip) , (25) 



T on (p) Air V a 2 

with 1/a -C Ao , and r e ~ 1 / Ao , where Ao is the "natural" scale set by the theory underlying the low- 
energy effective theory. Matching expressions ( p4[ ) and ( p5| ) for Ao > p > 1/a it is straightforward 
to deduce [||, || that the coefficients scale as 

Co ~ TT^ ° 2n ~ M^Af (26) 

provided that the scale /i is kept greater than, or of the order of, the scale 1/a. As pointed out by 
Kaplan et al. we may now choose the scale /i to be of order p. 

With this choice all loops involving the operator Co are of the same order and so the operator 
Co(V^V') 2 must be treated non-perturbatively. This power counting can be understood from the 
viewpoint of the renormalization group (3|, ^] as an expansion around a nontrivial IR fixed point. 
This fixed point describes systems with a bound state at exactly zero energy. In the PDS scheme, 
the potential at the fixed point is independent of energy or momentum and scales like 1/fJ,, as in 
Eq. @. 
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From Eq. pq ) we see that the effects of the operators C2WP are suppressed by a factor (p/Ao) n 
relative to the non-perturbative effects of the dimension-six operator. Therefore, these higher- 
dimensional operators may be treated perturbatively. Hence we can re-expand the first term in 
Eq. in the form 



1 



C2 2 , 

cf + 



C2- 2 

cl 



C4 
cl 



P 4 + 



(27) 



The effective range expansion (|25| ) can then be reproduced to any desired order in this p/Ao 
expansion. 

In an RG approach the coefficients in the effective range expansion are in one-to-one correspon- 
dence with the coefficients of the RG eigenfunctions in the expansion of the potential around the 
fixed point. These coefficients scale with definite powers of the renormalization scale fi in the PDS 
scheme or the cut-off in a Wilsonian approach |21j| . Although the scaling behavior is the same in 
the two approaches, the detailed form of the potential is not: the fixed-point and RG eigenfunctions 
all have a more complicated energy dependence when a cut-off is used. 



3.2 Scaling in PDS in a natural theory 



Note that in a natural theory, where 1/a ~ Aq, matching (25) to (pj) in fact yields 



G 



1 



2n 



MA 2 n+1 



(28) 



provided only that \x < Aq. Since the scale // does not enter these scaling relations we are free to 
choose (j, = 0, which is, of course, the choice of minimal subtraction. However /1 / with p ~ fj, 
is also a choice that leads to good power-counting, Clearly the tree-level effects of the operators 
Cm are suppressed by a power (p/Ao) 2n . Moreover, the effect of a loop with superficial degree 
of divergence 2n + 1 is down by a factor of (p/Ao) 2n+1 , relative to the tree-level Co- Thus in the 
natural theory there is no need to treat the effects of Co non-perturbatively, and the effects of 
higher-derivative operators are more strongly suppressed than in the unnatural case. 



3.3 Scaling in uvPDS 




Figure 3: Characteristic loop graph. 

We can now deduce the effect of modifying the definition of I n by adding the terms for m = 1 to 
m = n which appear in Eq. (^). Recall that if a cutoff regulator is used these terms correspond to 
cubic and higher power-law divergences. Attempting to solve the theory exactly, as was done above 
for the PDS subtraction, becomes much more complicated (see Refs. fi"5| , 16, ^] for examples), as 
each integral I n has a different number of subtractions. This complexity gives rise to nonlinear 
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relations between the bare and renormalized coefficients. However, the crucial point in any such 
attempt is that if we keep the scale [i well below the scale Ao then the effect of all power-law 
divergences will be suppressed by powers of (x/Aq, and so they make no change to the scaling of 
coefficients quoted above. For instance, if we examine a loop as shown in Figure | with one insertion 
of C 2 ( n _n and one of C21 then, considering for the present only the pieces of the sum in Eq. (p0[ ) 
corresponding to power-law divergences of degree three and above, the loop becomes 

M n ///\ 2m+1 1 

m=l 

Given the scaling above for the coefficients C2j the mth term here produces scaling of the coefficient 
C2(n-m) of the form 



Compared to the leading scaling of this coefficient C 2 ( n _ m ) ~ 1 / (M fi n ~ m+1 A^'" 1 ) the scaling ( j30|) 
is suppressed by a factor of (/x/Ao)" 1 . Thus, even if the terms in Eq. ( p0|) which correspond to 
power-law divergences of higher degree than those considered in PDS were included they would not 
modify the PDS scaling (|26|). 

By contrast, the terms arising from poles in D = 5 and above, which sum up to produce the 
logarithm of Eq. (p^) , would modify the PDS scaling (|26|) if they were included in the definition of 



I n . To see this let us only retain the terms of the sum in Eq. (20) for m < 0. Then 



2n M f : _ , ix fn + 2p 



so solving the Lippmann-Schwinger equation gives 



1 - 1 + «e_** tog f/i±5') + "*. ( 32) 



Expanding out the log in powers of p/fi it is straightforward to see that matching to the expression 



|25|) requires 

1 

Ml? 

This, of course, is the same scaling one would expect in a cutoff theory in which the cutoff was kept 



below the scale A [11, 17]. 



Now, if the coefficients do scale in this way then for momenta p ~ /i the effect of the higher- 
derivative operators C2 n p 2n is not suppressed relative to that of the "lowest-order" operator Cq. 
Thus if uvPDS is used to regulate the divergent loop integral the resulting power-counting for the 
C's is such that one cannot justify any truncation of the sum over all higher-derivative operators. 
Consequently, we conclude that the choice ( |20| ) does not lead to good power-counting for the 
coefficients C2 n - 

The physical reason this occurs is that the inclusion of terms that contain negative powers of 
// in the definition of the integral J" vPDS introduces energy-dependence in the amplitude at scale 
fi. Since /i is much less than the natural scale of energy dependence Ao, this energy dependence 
is unnaturally rapid. It must be cancelled by using larger values of the coefficients C2 n - In fact 



as shown in [21], this occurs automatically within an RG treatment where the unnatural energy 



9 



dependence is absorbed into the forms of the fixed-point potential and the RG eigenfunctions for 
the expansion around it. In our case we would have the fixed-point potential 



"M=|-^ + SW*±£]) ■ (34) 



4vr 2 4vr 2 \n~ 2p, 

If this potential is iterated via the Lippmann-Schwinger equation and uvPDS (which we have already 
shown is equivalent to a cutoff) applied to define the divergent integral then all the spurious energy- 
dependence introduced by the terms in the sum for I n with negative-powers of [i are cancelled. The 
analysis given here shows that this is essentially equivalent to just ignoring the effects of the poles 
in D > 5 which lead to logarithmic behavior in the inverse amplitude altogether. 

One difference between uvPDS (or a cut-off) and PDS a la KSW is that it distinguishes between 
energy and momentum dependence in the potential. The RG eigenfunctions that correspond to 
the terms in the effective range expansion are purely energy-dependent [21]. Ultimately they can 
be thought of a terms in an energy-dependent pseudopotential that acts as an energy-dependent 
boundary condition on the wave function at the origin §, pi. There are in addition momentum- 
dependent eigenfunctions with different scaling behaviors, but these do not contribute to the on-shell 
scattering amplitude. 



4 Summary 

The integral /„ is a key ingredient of non-relativistic effective field theories in general, and EFTs 
of the nucleon-nucleon interaction in particular. Since this integral has infinitely many poles in the 
complex space of dimensions there are infinitely many ways to define it, depending on how many 
of these poles we choose to subtract in our redefined I n . The original suggestion of Ref. || which 
removes the pole at D = 3 therefore exists amidst myriad alternative schemes. The alternative 
corresponding to removing all of the poles which arise due to ultraviolet divergences of I n exactly 
reproduces the result obtained by simply regulating I n via a cutoff. 

We have also analyzed the way that coefficients in the effective field theory Lagrangian scale in 
regularization schemes using these more general definitions of I n . This analysis shows that, in the 
presence of an unnaturally large scattering length, PDS is the simplest definition of a subtracted 
I n which leads to good power-counting for the coefficients that appear in the effective field theory 
Lagrangian. 
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A "Maximal" subtraction 

Here we consider the role played by the poles in D < 1 — 2n. As mentioned above, these correspond 
to infrared singularities of the original integrand. Since the integral of interest has no infrared 
divergences in D = 4, this discussion is not particularly relevant to nature. However, it leads to an 
intriguing result, and so we present it in this appendix. 
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To cancel the poles of T( 2n+ ^ 1 ) in D < 1 — 2n dimensions we must add a term 

oo / x 2m+l -i 



2vr 2 ^„\2j L>-3 + 2m' 

m=n+l x * 

We can then define a scheme which we call maximal subtraction, that subtracts all poles in the 
complex D-plane: 

M 00 /n\2m+l l 

m=— 00 v ' 

in all dimensions -D in the complex D-plane. We now specialize to the case D = 4. 

The sum in Eq. (^6[) is divergent for any value of the momentum p. In order to define a sum 
for this series we first define the function f(p) via 

M 00 / ;i\ 2m+l -1 
m=0 

If jLt/2 < p then this series may be summed to yield 

In the region of the complex p-plane, fj,/2 < \p\ < R, with R some large real number, this function 
is an analytic function of p. Hence we can define a function f(p) which is the analytic continuation 
of f(p) into the region < \p\ < R. We make the analytic continuation 



p + p/2 



2n iM P 



+ P ^. (39) 



p-p/2 

We now define to be the result found if the sum over positive m in Eq. (|36|) is replaced 



by its analytic continuation /, thus: 



fMaxS _ f f,.\ , J2n^ M P 



Taking the limit D — > 4 and evaluating the (finite) integral I n {p) implies 

J^ axS = 0. (41) 

By subtracting all the poles in the complex D-plane we have forced a definition of the value of the 
divergent integral I n which gives zero. 
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